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ABSTRACT 
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We  are  concerned  here  with  the  deflagration  of  reactants  that  are  produced  by 
gasification  at  the  surface  of  a solid  or  liquid.  In  particular,  the  influence  of 
pressure  on  the  gasification  rate  is  examined  under  various  conditions.  In  contrast 
to  earlier  treatments,  the  basic  phenomena,  including  several  unsuspected  ones,  are 


clearly  uncovered  by  activation-energy  asymptotics. 
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SIGNIFICANCE  AND  EXPLANATION 


The  burning  of  a solid  actually  consists  of  two  processes: 
gasification,  or  pyrolysis,  of  the  solid  at  its  surface  and  reaction 
of  the  gases  so  produced  to  form  a flame.  The  same  two  processes  are 
present  when  a liquid  hums,  the  gasif icatian  now  being  evaporation. 
Here  we  are  concerned  primarily  with  the  effect  of  pressure  on  the 
gasification  rate,  in  particular  conditions  under  which  the  flame  is 
extinguished  (or  ignited).  Radiation  from  (and  to)  the  surface  is 
important,  as  is  heat  exchange  between  the  gases  and  their  surroundings 
(by  radiation  or  conduction) . 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


MATHEMATICAL  THEORY  OF  LAMINAR  COMBUSTION 
IV.  STEADY  BURNING  OF  A LINEAR  CONDENSATE 

J.D.  Buckmaster  and  G.S.S.  Ludford 

I 

1.  Responses. 

In  Secs.  II. 3, U several  methods  of  supplying  a reactant  at  a station  x = 0 were  con- 
sidered.  In  each  case  four  relations  between  the  three  surface  values  (II. 2k),  M ana  D 
are  prescribed  and  then  the  flame  eigenvalue  determines  them.  For  the  switch-on  reaction 
the  surface  values  themselves  are  prescribed  along  with  D and  then  the  eigenvalue  gives  M. 
For  the  flameholder  M,  D,  Ts  and  Js  = 1 are  prescribed  and  the  eigenvalue  gives  Tfc, 
i.e.  Yg  + Ts.  (The  third  example,  namely  the  vaporizing  liquid,  is  discussed  at  the  end  of 
Sec.  7.) 

The  switch-on  reaction  provides  the  simplest  example  of  a response.  With  the  surface 

values  (i.e.  the  switch-on  temperature  Tg  and  the  upstream  state  Y ^.T  J)  held  fixed,  the 

flame  speed  (i.e.  M)  is  determined  as  a function  of  pressure  (i.e.  D).  The  response  curve 
-2 

is  a parabola  DM  = const.,  where  the  constant  is  determined  by  T = Y + T and  J = Y 

00  —00  —00  S —00 

Note  that,  while  x,  is  fixed,  the  location  of  the  flame  varies  along  the  parabola  because 
M is  involved  in  the  length  unit. 

For  the  flameholder,  where  Tg  and  Jg  are  fixed,  the  flame  temperature  T^  becomes 
a function  of  both  M and  D.  Such  a function  can  be  described  by  its  sections 
MT>)2exp(e/2T<>))  = vS)/0  = const.,  which  are  parabolas  flattened  along  the  T^-axis  near  the 
origin.  However,  not  all  of  the  curve  represents  the  response:  T^  must  lie  between  Tg 
and  Tg  +1  if  x„  is  to  lie  in  the  range  (0,®).  These  limits  correspond,  respectively,  to 
all  and  none  of  the  heat  released  at  the  flame  being  conducted  back  to  the  supply.  When 
T = T the  flame  has  reached  the  surface  and  the.  analysis  of  Sec.  II. 5 shows  that  M may 
then  be  decreased  indefinitely  without  changing  T^.  When  T^  = Tg+1  the  flame  has  become 
remote  and  the  same  section  shows  that  M may  then  be  increased  indefinitely.  The  resulting 
response  is  similar  to  that  in  Fig.  1 with  M,D  replaced  by  T^,M.  In  practice  the  remote 
flame  signals  an  effective  extinction  (cf.  Sec.  2). 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024 . 
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FIG.  1 


Solid  pyrolysis  under  adiabatic  conditions: 
•■••stable  response. 


basic  M,D-curve; 


A more  complicated  response,  again  of  M to  D,  occurs  when  gasification  at  the  surface 

of  a solid  or  liquid,  lying  in  x < 0,  supplies  the  reactant.  We  still  have  = 1 and 

T (=Y  + T ) can  be  calculated  from  an  overall  enthalpy  balance.  Under  adiabatic  conditions, 
oo  s s 

T«,  itself  is  determined  but  otherwise  it  is  expressed  in  terms  of  Tg  and  M.  A further 
relation  comes  from  specifying  the  nature  of  the  gasif4  .ion : in  pyrolysis,  the  surface 
temperature  Tg  determines  the  production  rate  M;  in  vaporization,  it  determines  the  partial 
pressure  represented  by  DYg . The  fourth  relation  is  prescription  of  D,  and  then  the  flame 
eigenvalue  is  an  equation  for  M,  the  T^  in  it  being  known  (implicitly)  as  a function  of 
M.  Our  task  is  to  determine  how  M varies  with  D when  the  remote  temperature  of  the 
condensate  and  the  background  temperature  are  held  fixed. 
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2.  Experimental  Results.  Extinction. 

If  a solid  or  liquid  contained  in  the  left  half  of  a long  straight  uniform  tube  gasifies 
into  a combustible  mixture  at  its  surface  then  a steady  state  may  be  attained  ir.  which  the 
surface  recedes  and  is  followed  at  a fixed  distance  by  a flame  separating  unburnt  gas  from 
burnt.  An  observer  moving  with  the  surface  sees  a constant  flux  M of  unburnt  gas  towards 
the  flame  and  an  equal  flux  of  burnt  gas  beyond.  The  problem  of  determining  the  combustion 
field  is  of  the  type  which  arose  in  Sec.  II. 3 while  resolving  the  cold-boundary  difficulty, 
if  the  combustible  mixture  can  be  treated  a single  reactant;  in  fact,  the  vaporizing  liquid 
treated  there  is  one  of  the  two  main  questions  here,  and  it  will  have  to  be  reconsidered. 

The  other  main  question  is  the  pyrolysing  solid. 

Since  the  aim  is  to  determine  the  response  of  M to  D,  it  might  be  thought  that 
equation  (11.22)  has  the  answer:  the  response  is  parabolic.  However,  the  condensate  has  two 
important  effects.  First,  it  fixes  the  position  of  the  flame,  which  must  be  somewhere  between 
its  surface  and  infinity.  The  use  of  the  formula  (11.22)  is  thereby  restricted  and  indeed  when 
the  flame  approaches  the  surface  or  recedes  to  infinity  it  must  be  replaced  by  the  appropriate 
formula  from  Sec.  II. 5.  Secondly,  it  can  make  T dependent  on  M or  D and  hence  change 
completely  the  response  predicted  by  the  original  formula.  These  two  effects  are  the  subject 
of  the  present  chapter. 

Experiments  which  approximate  the  situation  described  above  have  been  done  for  solids 
which  gasify  by  surface  pyrolysis  (Johnson  & Nachbar , 1962;  Nir , 1973).  Three  basic  phenomena 
are  found:  (i)  the  burning  rate  increases  with  ambient  pressure,  but  very  little  beyond  a 
certain  pressure;  (ii)  there  is  a minimum  pressure  beyond  which  gaseous  deflagration  will 
not  take  place  (extinction);  and  (iii)  preheating  the  solid  lowers  this  minimum.  All  three 
effects  are  exhibited  by  the  theory,  even  when  absorption  or  emission  of  radiant  energy  by 
the  surface  and  heat  loss  or  gain  through  the  wall  of  the  tube  are  neglected.  These  depar- 
tures from  adiabatic  deflagration  will  be  considered  later;  here  we  Just  note  the  partial 
analysis  given  by  Johnson  & Nachbar  (1962)  based  on  distributed  heat  loss  from  the  condensed 
phase . 

The  analysis  follows  Buckmaster,  Kapila  h Ludford  (1976)  and  reveals  that  extinction, 
by  which  is  meant  the  disappearance  of  gaseous  deflagration  and  not  the  cessation  of  pyrolysis, 
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can  be  of  two  kinds.  In  true  extinction,  the  flame  suddenly  disappears  as  the  pressure  is 
lowered.  However,  the  pyrolysis  continues  at  a lower  temperature,  the  products  passing  away 
unburnt  in  what  is  known  as  flameless  combustion.  Effective  extinction,  followed  by  flame- 
less combustion  without  sudden  drop  in  surface  temperature,  occurs  when  the  flame  recedes 
to  infinity  so  that  the  products  of  continuing  pyrolysis  effectively  pass  away  unburnt.  The 
disappearance  of  a flame  is  generally  taken  to  be  true  extinction,  although  Nir  comes  close 
to  recognizing  effective  extinction  in  his  experiments. 

We  shall  find  that  extinction  is  always  of  the  effective  kind  for  adiabatic  deflagration, 
and  may  be  so  for  radiative  loss  or  gain  (when  a separate  phenomenon  of  ignition  is  possible). 
We  shall  also  show  that  effective  extinction  is  changed  to  true  extinction  by  heat  losses 
through  the  tube  wall  (Kapila  & Ludford,  1977). 

Apparently  no  experiments  have  been  reported  on  vaporizing  liquids.  A more  realistic 
treatment  than  that  presented  in  Sec.  II.  1*  reveals  the  striking  effects  detailed  later,  which 
deserve  experimental  verification. 

3-  Solid  Pyrolysis  Under  Adiabatic  Conditions. 

The  gas  mixture  will  be  treated  as  a single  reactant  with  Lewis  number  1.  The  object  is 
to  determine  the  combustion  field  in  general,  and  M in  particular,  given  the  remote  tem- 
perature T M of  the  solid  and  the  pressure  level  in  the  gas  phase,  i.e.  D.  For  that 
purpose  (see  Sec.  II. 4)  we  need  only  determine  (in  addition  to  M)  the  values  Ys,Tg 

» I 

and  Yg  = -Tg  which  the  solid  in  x < 0 presents  to  the  gas  phase  in  x > 0. 

_2 

In  fact,  to  determine  the  eigenvalue  A=  H D it  is  only  necessary  to  know  the  mass 

f 

flux  fraction  J = Y - Y and  the  flame  temperature  T = T_  + Y . If  none  of  the  products 
sss  »ss 

of  the  gaseous  reaction  is  produced  or  absorbed  at  the  surface  of  the  solid,  then  = 1 
and 

(1)  M"2D  = 82exp(e/Tj/2T^  , 

a relation  determining  M once  T—  is  known.  The  formula  has  a mixed  character:  6 and  T 

n 

have  no  dimensions  while  M and  D have  those  of  mass  flux.  Several  such  equations  will 
appear  in  the  sequel.  Note  that  the  presence  of  the  solid  (or  whatever  else  is  creating  the 
stream  of  reactant)  is  felt  solely  through  T^  , a fact  first  recognized  by  Williams  (1973). 
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T = Y + T can  be  obtained  without  calculating  surface  values.  In  the  absence  of 
co  s s 

radiation  to  or  from  the  surface  and  of  losses  elsewhere,  the  flux  of  enthalpy  at  x = + °° 
must  equal  that  et  x = (since  the  kinetic  energy  is  negligible).  Reverting  to  dimensional 
quantities,  we  may  write 

(2)  (h°  - Q)  + c (T  - T°)  = h°  + c (T  - T°), 

p 00  P —00 

where  h°  - Q is  the  heat  formation  of  the  product  of  the  deflagration  and  A 

denotes  values  for  the  solid.  Solving  for  T and  writing  the  result  in  dimensionless  form 

yields 

(3)  T^  = + (1  -k)T°  + q + 1, 

where 

(U)  k = c /c  and  q = (h°  - h°)/Q. 

P P 

Here  q,  the  heat  of  pyrolysis,  is  the  difference  between  the  heat  of  formation  of  the  solid 

and  that  of  the  gaseous  reactant,  both  at  the  standard  temperature  T°,  expressed  in  units  of 

the  heat  of  combustion  Q.  Equation  (3)  shows  that  T (which  we  shall  suppose  to  be  positive) 

is  independent  of  M and  is  a linear  function  of  the  remote  temperature  of  the  solid. 

Fig.  1 sketches  the  curves  (l)  for  two  values  of  T m. 

A third  condition  must  be  added  to  J = 1 and  T + Y = T in  order  to  determine 

s s s °° 

all  the  surface  values;  that  takes  the  form  of  a pyrolysis  law.  The  rate  of  gasification  is 
supposed  to  depend  only  on  the  surface  temperature  and  not,  for  example,  the  pressure.  We 
shall  take 

(5)  M = kTs  exp(-0/Tg), 

where  k and  0 are  constants.  Note  that  the  pyrolysis  occurs  at  any  surface  temperature, 
so  that  the  term  extinction  (or  ignition)  must  refer  to  the  gaseous  deflagration,  as  we  have 
supposed.  The  law  should  be  viewed  as  determining  the  surface  temperature  T,.  which  produces 
the  required  flux  M into  the  gase  phase.  Various  modifications  have  been  suggested,  all  of 
which  lead  to  similar  results.  The  only  essential  feature  is  that  M should  be  a steadily 
increasing  function  of  Tg.  Apart  from  its  simplicity,  the  law  adopted  here  has  the  advan- 
tage of  firm  theoretical  foundations  in  gas  kinetics  and  reactive  solids.  Its  role  in 


f 
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Fig.  1 is  just  to  make  Tg  a parameter  on  the  curves. 

The  gas  phase  is  now  completely  determined,  without  reference  to  the  conductivitiy  of 

the  solid.  Only  the  enthalpy  of  the  solid  enters  into  the  result  (l);  and  the  remaining 

details  of  the  gas  phase  come  from  adding  the  pyrolysis  law.  The  role  of  the  conductivity 

is  to  fix  the  distribution  of  temperature  between  T and  T so  that  the  total  heat  flux 

-oo  s 

at  each  point  of  the  moving  solid  is  constant. 

Whatever  the  value  of  8 , the  surface  temperature  must  lie  within  certain  bounds; 

Too  = Ys  + Ts  and  0 < Yg  < 1 imply 

(6)  T - 1 < T < T . 

00  s 00 

For  finite  8,  it  is  not  immediately  clear  that  a solution  is  then  guaranteed;  but  in  the 
limit  0-*»  such  is  the  case.  Only  the  location  (II. lU)  of  the  flame  need  be  checked.  If  we 
note  that 

(7)  Ts  = -Ys  = 1 - Ys  = 1 + Ts  - T-» 

a result  expressing  heat-flux  balance  in  the  gas  phase,  then 

(8)  x#  = -*p(l  + Tg  - Tj 

and  the  inequalities  (6)  ensure  0 < x#<  «°  ,as  supposed.  It  is  a characteristic  ieeture  of 

activation-energy  asymptotics  that  such  questions  can  be  settled  simply. 

Fig.  1 shows  the  corresponding  end  points  Pg  and  , so  designated  because  at 

P the  flame  has  moved  back  to  the  surface  and  at  P has  receded  to  infinity.  These 

s 00 

two  possibilities  were  dealt  with  before  (Sec.  II. 5),  requiring  reconsiderations  of  the 
structure.  It  was  found  that,  with  conditions  at  the  surface  fixed,  the  remote  flame  allows 
D to  be  decreased  to  zero  while  the  surface  flame  allows  it  to  be  increased  indefinitely, 
both  without  change  in  M;  that  is  just  what  the  pyrolysis  law  demands.  Thus  the  M,D-curve 
is  completed  by  horizontal  lines,  one  stretching  from  Pg  to  D = ■»  and  the  other  from  P^ 
to  D = 0.  Of  course,  if  T is  less  than  1,  the  left  inequality  (6)  is  ineffective,  i.e. 

P does  not  exist  and  the  curve  extends  to  the  origin.  The  following  remarks  will  apply 

when  T^  > 1 and  must  be  slighly  modified  otherwise. 
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The  response  curve  now  shows  that  (i)  M increases  with  D,  but  not  beyond  ; and 
that  (ii)  there  is  a minimum  D,  corresponding  to  P , before  which  there  is  effective 
extinction.  The  third  phenomenon  (iii)  listed  in  Sec.  2 concerns  the  change  ir.  D,  for 


T = - 1,  as  increases  linearly  with 


exp(-6/Tj  < 0 


We  find 


with  relative  error  0(0  ),  since  dM/dT  = dM/dT  is  0(l).  The  effect  (iii)  is  clearly 

<»  s 

present  and  is  caused  by  the  exponential  factor,  which  decreases  as  T increases  irrespec- 
tive of  how  M changes. 

h.  Radiation  from  the  Surface. 

It  is  of  interest  to  see  how  the  picture  is  modified  by  radiation  between  the  surface  of 
the  solid  and  the  surroundings.  Distributed  exchange  will  be  treated  later.  If  the  background 
temperature  is  , there  is  a heat  loss  (or  gain  if  negative) 

(10)  r=o(T^-Tj), 

where  o is  a positive  constant  which  is  not  always  small  in  practice.  Only  Spalding  (i960) 
has  treated  radiative  exchange  with  surroundings  at  nonzero  temperature  (T^  9*  0),  and  then 
in  an  approximate  fashion  which  precluded  heat  gain  (T^  > Ts).  As  a consequence 

he  missed  the  most  striking  results.  Experiments  on  such  loss  have  been  reported  by  Levy 

& Friedman  (1962). 

The  loss  changes  the  overall  enthalpy  balance  to 

(11)  T^  = T*  - r/M> 

where  T^  is  the  flame  temperature  under  adiabatic  conditions,  given  by  the  formula  (3); 
but  nothing  else  changes.  T is  now  a function  of  M,  both  directly  and  through  T in 

co  S 

r;  inserted  in  the  response  formula  (l),  it  provides  a D which  is  no  longer  proportional 

to  V? . Note  that  the  pyrolysis  law  now  affects  the  shape  of  the  M,D-curve. 

Consider  first  T^  = 0.  To  determine  the  new  response  curve  it  is  helpful  to  sketch 
the  graph  of  T versus  T for  fixed  T , and  that  has  been  done  in  Fig.  2.  There  is  a 

* co  S -00 

maximum  at  E,  corresponding  to  T = 9/3,  from  which  the  descent  is  monotonic  to  -«»  as 


V 


Ts  -*■  0 or  oo.  The  shape  of  the  M,D-curve  is  obtained  by  rotation  through  90°  since  the 
asymptotically  correct  result 

(12)  D-1dD/dM  = ~(6dT  /T2dM)dT  /dT 

S oo  oo  s 

shows  that  all  slopes  have  their  signs  changed.  Fig.  3 gives  sketches  of  the  resulting 
C-shaped  curve,  the  leftmost  point  corresponding  to  E. 

As  in  the  absence  of  radiation,  not  all  of  the  curve  is  acceptable.  In  Fig.  2 only  the 
portions  lying  between  the  lines  = 1_  and  T = Tg  + 1 satisfy  the  inequalities  (6), 

and  these  depend  on  T . (An  increase  in  T merely  translates  the  curve  upwards.)  If 
T is  too  small,  there  is  no  such  portions;  otherwise  five  possibilities  arise  depending 
on  whether  the  curve  intersects  one  or  both  of  the  lines  with  E between  them  or  not.  These 
lead  to  the  five  parts  of  Fig.  3,  where  progress  is  from  (a)  to  ( bor  b')  to  (c)  to  (d)  as 
T increases.  In  each  case  the  end  points  of  the  portion  are  marked  to  correspond  with 
the  position  of  the  flame,  and  horizontal  completions  are  included.  The  complete 
responses  are  shown  by  unbroken  lines. 
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It  ir  widely  believed  that  if  M decreases  as  D increases  the  combustion  is  unstable 
(Emmons,  1971)  and,  since  there  is  some  analytical  evidence  to  support  the  idea,  we  shall 
adopt  it  here.  (Surface  burning  corresponding  to  a horizontal  extension  from  such  a point 
will  also  be  assumed  unstable).  The  whole  of  Figs.  3a  and  V are  then  discarded,  the  pre- 
heating being  too  weak  for  stable  combustion.  The  surviving  parts  of  Figs.  3b,  c and  q 
are  drawn  heavily,  as  the  responses  to  be  expected  in  an  experiment.  Clearly  they  exhibit 
the  phenomena  (i)  and  (ii)  mentioned  in  Sec.  2. 

The  minimum  T ^ for  stable  burning  can  be  found  by  setting  T^  = Tg  = 0/3  in  (ll) 
and  solving  for  T As  the  preheating  increases,  the  true  extinction  of  Figs.  3b  and  c 

will  occur  until  T reaches  the  value  given  by  (ll)  for  T - 1 = T = 6/3.  For  all 

—oo  oo  S 

larger  values  of  T ^ the  effective  extinction  of  Fig.  3d  will  occur,  as  under  adiabatic 
conditions.  The  extinction  value  of  D is  given  explicitly  by  the  formuLa  (:  , where  M 
and  T^  are  to  be  calculated  from  (5)  and  (ll)  by  setting  Tg  = 9/3  for  true  extinction 
and  T = T - 1 for  effective  extinction.  The  third  effect,  that  the  extinction  value 

S oo 

of  D decreases  as  T increases,  then  follows  from  the  derivative  (9)  since  T still 
increases  with  T 

—00 

$■  Background  Radiation. 

If  T,  / 0 in  the  radiation  term  (10)  several  new  responses  can  occur.  These  are  not 
b 

associated  with  Tg  large  (where  the  results  of  Sec.  1 apply)  but  with  Tg  small,  which 

not  only  allows  to  change  the  sign  of  r but  also  enhances  the  latter's  effect  in  (ll) 

through  the  smallness  of  M.  The  left-hand  side  of  the  curve  in  Fig.  2 is  thereby  bent 

upwards  to  give  Fig.  4.  Whether  the  result  is  a monotonic  curve  (Fig.  la)  or  one  with  a 

minimum  and  a maximum  Eg  (Fig.  lb)  depends  on  the  size  of  T^:  the  polynomial 

factor  3T^  - 9T1*  + T^T  + 9T.!*  in  dT  /dT  has  lust  two  positive  zeros  (both  lying  between 
S S D S D 00  S 

and  9/3)  if  is  less  than  0.l68§  (<§/3)  and  no  positive  zero  otherwise.  Since  the 

curve  must  cross  the  strip  (6)  there  is  always  deflagration  (albeit  unstable)  for  some 
range  of  D,  in  contrast  to  Sec.  U where  T ^ had  to  be  large  enough. 

The  monotonic  curve  for  T^  > 0.1689  leads  to  a response  which  is  qualitatively 
the  same  as  for  adiabatic  conditions  (Fig.  l).  The  strong  background  radiation  tends  to 
compensate  for  losses  from  the  surface. 
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The  minimum-maximum  curve  for  <0.1666  is  cut  by  l5'-lines  in  at  most  three 

points.  Such  curves  present  30  possibilities  depending  on  how  the  lines  = To , T + 1 

are  cut (times  and  order)  and  on  whether  E , E0  lie  above,  between  or  below  them.  Three 
of  these  can  be  ruled  out  because  the  present  curve  has  only  one  inflexion  point;  others 
lead  either  to  unacceptable  responses  or  to  responses  which  have  already  appeared  for  = C . 
The  remainder  fall  into  five  groups,  according  to  similarity  of  response,  and  one  member 
(the  simplest)  from  each  of  the  five  will  be  presented  here.  Fig.  5 shows  these  responses, 
i.e.  the  M,  D-curves  with  parts  deleted  by  the  limitations  (6)  and  the  stability  require- 
ment of  positive  slope  but  completed  by  horizontal  lines.  The  first  four  of  these  responses 
have  two  branches,  each  ending  on  the  left  with  either  true  or  effective  extinction.  In 
(a)  and  (b)  both  branches  have  a plateau  to  the  right  and  we  shall  see  that  the  plateau  on 
the  upper  branch  in  each  of  (c)  and  (d)  is  reached  for  D large  enough.  The  three  phenomena 
in  Sec.  2 are  therefore  exhibited  once  more.  The  sole  exception  is  (e),  which  has  no  upper 
branch. 

We  shall  now  construct  the  probable  sequence  of  events  first  as  D increases  from  0 
to  00  and  then  as  D decreases  from  °°  to  0.  When  there  is  only  one  effective  extinction 
on  the  left,  as  in  (a),  (c)  and  (e),  the  lower  branch  is  followed  as  D increases.  The 
same  is  true  for  (b)  and  (d)  if,  as  we  shall  assume,  weaker  burning  is  preferable  to  strong 
when  the  latter  is  not  already  established.  In  (a)  and  (b)  the  lower  branch  will  be  followed 
all  the  way  through  surface  burning  to  infinity  but  in  (c)  and  (d)  there  will  be  a JumF  to 
the  upper  branch,  i.e.  ignition,  as  the  rightmost  point  is  approached  on  the  lower.  As  D 
decreases  in  (c)  and  (d)  the  upper  branch  will  be  followed,  for  (d)  all  the  way  to  effective 
extinction  but  for  (c)  with  a Jump  to  the  lower  branch  as  true  extinction  is  approached  on  the 
upper.  In  (a)  and  (b)  we  invoke  the  weaker  burning  assumption  to  see  that  the  lower  branch 
will  be  followed  all  the  way  to  effective  extinction.  In  (e)  there  is  true  extinction  as 
D increases  and  ignition  as  it  decreases.  After  extinction  and  before  ignition  there  is 
no  steady  burning  if  we  persist  in  believing  that  the  surface  burning  which  emanates  from 
the  unstable  branch  is  itself  unstable. 

The  above  description  is  based  on  several  unproved  p.  ’ nciples:  no  part  of  an 
M,D-curve  with  negative  slope  can  be  attained  physically  because  the  combustion  is  unstable; 
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a branch  will  be  followed  as  far  as  it  can  be  as  D increases  or  decreases;  and  when  a 
choice  remains  weaker  burning  will  occur.  Whether  the  other  branches  in  (a)  and  (b)  can 
ever  be  attained  is  an  interested  open  question  but,  supposing  they  cannot,  (c),  (d)  and 
(e)  are  the  only  really  new  responses.  Even  then,  only  (e)  fails  to  exhibit  the  three 
basic  phenomena. 

However,  data  of  Johnson  & Nachbar  (1962)  and  Guirao  & Williams  (1971)  for  ammonium 
perchlorate  can  be  interpreted  as  lying  on  the  strong  burning  branch  of  Fig.  5(a) , see 
Buckmaster,  Kapila  & Ludford  (1976),  suggesting  that  the  weak  burning  branch  is  not 
applicable.  That  would  be  the  case  if  pyrolysis  ceased  below  a certain  temperature,  so 
that  the  law  (5)  only  applied  above  and  M was  zero  below.  The  whole  lower  branch  could 
< hen  be  eliminated,  leaving  only  the  upper  branch  for  the  combustion  to  follow. 

It  should  be  noted  that  the  bending  back  of  the  M,  D-curve  from  a C-shape  into  an 
S occurs  for  any  non-zero  T^  < 0.1689.  For  T^  > 0.1689,  the  S is  pulled  cut  even 
further  into  a monotonic  curve.  Even  weak  background  radiation  can  therefore  result  in 
quite  different  responses.  Far  from  being  negligible,  it  can  be  the  dominant  effect  under 
suitable  conditions.  The  apparent  contradiction  as  T^  -*•  0 is  due  to  the  nonuniformity 
of  the  limit.  However  small  T^  is,  there  a re  smaller  values  of  Ts , for  which  background 
radiation  changes  the  heat  loss  into  a heat  gain.  Moreover,  the  effect  of  the  gain  on  the 
enthalpy  of  the  reactant  is  magnified  by  the  smallness  of  the  mass  flux  at  such  surface 
temperatures.  However,  it  is  unlikely  that  such  a limiting  behavior  could  be  observed, 
because  of  the  pyrolysis  cut-off  at  low  temperatures  mentioned  above. 

In  practice,  the  most  important  features  of  the  response  are  the  ignition  and  extinction 
values  of  D.  In  any  particular  case,  these  are  easily  determined  by  finding  the  points 
E. , E„  or  the  appropriate  intersection  of  the  T , T -curve  with  T = T + 1. 

Jl  C-  S 00  00  S 

6.  True  Nature  of  Effective  Extinction. 

The  object  of  the  present  section  is  to  determine  how  the  previous  results  are 
modified  when  0(1/9)  heat  is  lost  or  gained  (by  radiation  or  lateral  conduction)  through- 
out the  gas  phase.  For  simplicity  we  shall  start  by  supposing  that  the  solid  phase  is 
perfectly  insulated.  The  main  conclusion  is  that  all  effective  extinctions  are  changed 
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into  true  extinctions  bj  heat  lore,  the  explosive  regime  only  being  maintainable  as  the 
limit  for  vanishingly  small  heat  gain. 

For  each  point  of  the  relevant  response  curve  determined  above,  we  calculate  the  new 
D (due  to  heat  exchange)  for  the  j'jne  value  of  M.  The  general  effect  of  heat  loss  can 
be  seen  from  the  formula  (1).  If  the  flame  temperature  is  reduced,  M can  only  be  main- 
tained by  increasing  D:  in  order  to  maintain  the  rate  of  reaction  at  a lower  temperature 
the  concentration  of  reactant  must  be  increased,  i.e.  the  pressure  raised.  The  pyrolysis 
law  (5)  ensures  that  T£  is  unaffected,  from  which  it  follows  that  the  temperature  distribu- 
tion in  the  solid  and  T'  are  also.  Yg  and  do  change,  though  keeping  Ys  - \'s 

equal  to  1. 

The  problem  then  is  to  calculate  the  perturbation  of  the  flame  temperature,  which  Kapila 
& Ludford  (1977)  did  by  straightforward  matching.  However,  we  shall  follow  Sec.  III. 3 in 
calculating  it  directly  from  the  change  in  enthalpy  from  the  surface  up  to  (and  including) 
the  flame  sheet.  Rather  than  modifying  equation  (III. 28),  we  shall  derive  the  result 
ab  initio  for  the  present  circumstances. 

The  governing  equations  are 

(13)  £(Y,1)  = - fi(T,l)  + 6M'2$(T)  = exp (-9 /T ) , 

where  <p  is  the  function  (III. 30)  with  T^  replaced  by  T^.  Integration  of  the  first 
equality  between  0 and  x*  + 0 immediately  yields  the  flame-temperature  perturbation 

(lU)  t-  = ?f|x,tC-M"2  A(T)dX’ 

0 

since  the  surface  values  T , T'  and  Y - Y'  are  not  perturbed.  Both  terms  on  the 

s s s s 

right-hand  side  are  to  be  evaluated  to  leading  order  and,  for  that  purpose,  the  approximations 

(Ts  + T^(eX-l)  for  0 < x < x„, 

T„  + M~^$(Tj>)(x#  - x)]  for  x > x, 

are  used,  to  obtain 

x# 

(l6)  t = - «f2  where  4>  =*(T  ) + f <HT  + T’(eX-l))dx, 

00  CD  J S S 
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a result  that  should  be  compared  with  the  formula  (III. 35).  The  two  terms  ir.  { represen' 
heat  exchanges  of  the  unburnt  mixture  (between  the  surface  and  the  flame)  with  the  burnt 
mixture  and  surroundings,  respectively.  (The  burnt  mixture  later  exchanges  the  same  heat 
with  its  surroundings).  It  follows  that 

(17)  D = D0(M)e*M 

where  D ( M ) is  given  by  the  eigenvalue  (l). 
o 

The  last  three  sections  have  been  an  investigation  of  Dq(M),  which  involves  the 

dependence  (ll)  of  T^  on  M both  directly  and  through  the  pyrolysis  law.  The  exponent 
_2 

$M  also  depends  directly  and  indirectly  on  M,  with  the  additional  involvement  of  the 
function  4>  . It  is  therefore  a complicated  matter  to  describe  the  dependence  of  D on  M 
in  detail,  especially  when  $ is  left  arbitrary.  However,  if  4>  is  not  too  large,  the 
general  shape  of  the  response  will  be  maintained  except  near  where  x„  -*■  » in  t. 

What  happens  there  depends  on  whether  the  surface  is  hotter  or  cooler  than  the  background. 
When,  as  is  usually  the  case, 

(18)  T (M  ) > T. 

s OD  D 

the  temperature  in  the  gas  phase  for  M above  is  everywhere  higher  than  T^, 

so  that  is  positive.  Consequently  t is  positive  and  tends  to  +«°  as  M -*•  M 
because  x„  does.  The  response  lies  to  the  right  of  that  in  Fig.  1 , 3 or  5 and  bends 
round  as  M -*■  to  form  a C.  We  conclude  that  heat  loss  in  the  gas  phase  changes  an 
effective  extinction  into  a true  extinction. 

In  the  less  likely  event  that 

(19)  TS(MJ  < Tfc, 

the  temperature  in  the  gas  phase  near  the  surface  is  lower  than  T^  for  at  least  a stretch 
of  values  above  . It  is  easily  seen  that,  as  M -*•  M^,  the  resulting  heat  gain 

eventually  overwhelms  any  heat  losses  further  away,  so  that  t becomes  negative  (if  it  is 
not  already  so)  and  tends  to  -«> . The  response  therefore  lies  ultimately  to  the  left  of 
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that  in  Fig.  1 , 3 or  5 and  ends  at  D = 0.  Even  then  the  effective  extinction  may  te  tree*  ■ 

by  a true  extinction:  if  T (M  ) is  sufficiently  close  to  T.  , the  function  4 vi  11 
" s 00  C 

become  negative  vury  close  to  and  before  that  will  be  increasing  rabidly  (because  of 

x*),  so  that  an  S is  formed.  For  smaller  values  of  Ts<fO  the  S straightens  out  into 
a monotonic  response.  Thus  the  explosive  regime,  corresponding  to  the  horizontal  line 
through  P^,  is  seen  to  be  a limit  for  vanishingly  small  heat  gain. 

Kapila  & Ludford  (1977)  have  given  details  of  the  above  picture  for  the  linear  lav 
(III. 30a)  when  the  ambient  is  at  the  same  temperature  as  the  remote  solid  (Tfc  = T_^)  and 
there  is  no  radiative  exchange  between  it  and  the  surface.  The  conditions  (10)  and  (19) 
then  hold  for  exothermic  and  endothermic  pyrolysis,  respectively.  They  also  consider  less 
than  perfect  insulation  of  the  solid  and  find  the  termination  points  are  unaltered  though 
the  curves  are  somewhat  distorted.  (T.  is  necessarily  the  same  as  T now,  since  other- 

D — < oo 

wise  the  solid  exchanges  an  infinite  amount  of  heat.)  The  reason  is  that  heat  exchange  in 
th'-  ;clid,  which  adds  a term  to  ♦ , is  bounded  as  M •*  K and  therefore  cannot  affect  the 
unboundedness  of  $ . 

Distributed  heat  loss  from  the  solid  phase  alone  was  considered  by  Johnson  & Nachbar 
(1962)  as  a way  of  modifying  the  C-shaped  responses  they  obtained  analytically  when  there  if 
radiative  loss  from  the  surface.  Under  their  assumptions  effective  extinction  cannot  occur. 

7.  The  Vaporizing  Liquid  Under  Adiabatic  Conditions. 

With  one  important  exception,  the  analysis  follows  that  in  Sec.  3,  the  superscript  * 
now  referring  to  the  liquid.  (It  applies  equally  well  to  a sublimating  solid.)  However, 
since 

(10)  L = -q  + (1  -k)(T8  - T°) 

is  the  latent  heat  of  vaporization,  which  must  be  increasingly  positive  at  all  temperatures 
Tg  of  interest,  the  heat  q must  be  negative  and  the  ratio  k <1.  We  are  dealing  with  a 
strictly  endothermic  process  in  contrast  to  pyrolysis,  which  is  usually  exothermic.  The 
exception  .Just  mentioned  is  the  use  of  the  pyrolysis  law  (5),  which  must  be  replaced  by  the 
Clausius-Clapeyron  relation 
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(21) 


(6  > 0) 


Vc 


= kTs  exp(-e/Ts) 


in  calculating  boundary  values.  Here  8 =Y(l  — k )/(r—  l),  where  < has  the  definition  (Ua) 
and  Y is  the  ratio  of  specific  heats  in  the  gas.  The  new  law  relates  the  partial  pressure 
of  the  vapor  at  the  surface  to  the  temperature  there  and  is  applicable  below  the  critical  point 
whenever  the  vaporization  is  rapid  enough  for  thermodynamic  equilibrium  to  be  achieved.  The 
critical  pressure  is  invariably  very  large  so  that  we  shall  allow  pc  to  become  indefinitely 
large,  as  is  needed  for  the  Damkohler  numbers  involved  in  the  asymptotic  theory.  Equation 
(21)  should  be  viewed  as  determining  Tg  for  given  D (i.e.  pg)  so  that  once  more  the 
surface  temperature  is  a parameter  on  the  parabola  (l),  see  Fig.  6.  For  the  pyrolyzing  solid 
the  surface  values  depended  on  M;  here  they  depend  on  D. 

Suppose  T >1  (otherwise  the  discussion  is  modified  in  an  obvious  way).  The  point  P^ 
corresponding  to  the  left  end  of  the  range  (6)  then  exists  and  is  marked  in  Fig.  6;  the  point 
P , corresponding  to  the  right  end,  lies  at  infinity  (according  to  equation  (21))  because  Yg 
vanishes  when  Tg  = T^.  So  long  as  Tg  is  away  from  the  ends  of  its  range  the  response  is 
identical  to  that  for  a pyrolyzing  solid.  Near  the  ends,  however,  the  two  responses  differ 
markedly,  due  entirely  to  the  difference  between  the  two  laws  (5)  and  (21). 

As  T -►  T^  - 1,  the  flame  sheet  recedes  to  infinity  and,  since  Y -»■  1,  the  pressure 

p (and  hence  D)  tends  to  a finite  non-zero  value.  The  analysis  of  remote  flames,  given 
c 

in  Sec.  II. 5,  shows  that  M can  then  increase  indefinitely,  with  conditions  at  the  surface 
fixed,  without  change  in  D;  and  that  is  just  what  the  Clausius-Clapeyron  law  demands. 

The  liquid,  which  is  at  its  saturation  temperature  for  the  ambient  pressure,  must  evaporate 
completely  before  the  pressure  can  be  lowered.  Combustion  plays  no  role  because  heat 
transfer  from  the  gas  phase  to  the  liquid  has  ceased;  the  entire  heat  needed  for  evaporation 
is  supplied  from  the  remote  end  of  the  liquid  which,  since  the  evaporation  is  endothermic 
(L  > 0),  is  now  hotter  than  the  surface  (cf.  equations  (3),  with  T^  = Tg  + 1,  and  (20)). 

The  process  is  represented  in  Fig.  6 by  the  vertical  line  through  P^;  before  it  occurs 
there  will,  of  course,  be  effective  extinction. 

As  T -►  T the  flame  moves  to  the  surface  and,  since  Y tends  to  zero  like  (T  - T ), 

the  pressure  p (and  hence  D)  tends  to  infinity  like  (T  - T )-1.  When  the  temperature 
c s 
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difference  becomes  of  order  4 the  eigenvalue  (l)  has  to  be  replaced,  the  surface- flame  r 

—2  -2 

(11.27)  showing  that  M D now  behaves  like  (T  - T ) ' . Hence  M behaves  like  (T  - ' 

00  S oc  £ 

or  E ',  and  we  conclude  that  the  M,D-curve  bends  down  from  the  parabola  to  asymptote  the 
D-axis.  However,  the  later  portion  of  the  curve  does  not  provide  an  acceptable  response  .un- 
its slope  is  negative. 

As  for  the  pyrolyzing  solid,  preheating  the  liquid  (and  therefore  increasing  T ) 

00 

lowers  the  rressure  needed  to  sustain  a given  burning  rate.  The  effect  of  preheating  on 

follows  from  the  relation  (21)  on  setting  Y « 1 and  Tg  = T - 1 and  is  the  opposite 
of  that  for  solid  pyrolysis:  the  pressure  p^  (and  hence  D)  increases  with  T and  hence 
T . A hotter  liquid  requires  a higher  pressure  to  produce  the  sudden  complete  vaporization. 

The  vaporizing  liquid  was  used  as  an  example  in  Sec.  II. 3 and  the  treatment  there  can 
now  be  seen  as  an  incomplete  answer  to  a somewhat  different  question.  There  the  temperature 
was  supposed  to  be  uniform  throughout  the  liquid  phase,  i.e.  the  remote  temperature  was 
maintained  at  the  surface  value.  As  a consequence  T , instead  of  remaining  constant,  varied 
with  the  surface  temperaturet  T = T in  the  determination  (3)).  The  M,D-relation  (l), 

— no  S 

which  was  not  treated  thoroughly  there,  is  no  longer  parabolic  since  T^,  now  depends  on 
D through  the  Clausius-Clapeyron  relation.  (We  shall  not  go  into  details  since  they  are  of 
limited  interest.)  Finally,  the  formula  (20)  shows  that  the  requirement  (6)  is  simply 
0 < L < 1,  as  found  before. 

8.  Radiative  Exchange  at  Surface.  Distributed  Heat  Exchange. 

In  considering  radiative  exchange  at  the  surface  it  is  convenient,  as  in  Secs.k  and  5, 

to  sketch  the  graph  of  T versus  T for  fixed  T which  now  results  from  eliminating 

00  s 

M from  the  overall  heat  balance  by  means  of  the  M,D-relation  (l)  after  having  eliminated  I 

from  that  relation  by  use  of  the  equilibrium  condition  (21). 

Fig.  7 shows  the  T^T^-curve  for  the  three  possibilities.  (In  Figs.  7(a)  and  (b)  we 

have  assumed  Ta  e 1 but  the  curves  are  similar  for  Ta  > 1.)  In  all  cases  the  curve 
00  00 

originates  at  (Ta  , Ta)  and  asymptotes  the  T -axis,  passing  through  (T  ,Ta  ) when  that 

oo  00  oo  D co 

point  lies  above  the  line  = T . In  case  (b)  there  is  always  an  extremum  E in  the 
strip  (6),  but  in  the  other  two  cases  it  may  lie  outside. 
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■curves  for  vaporizing  liquid  with  radiation  . 


The  asymptotic  formula  (12)  still  shows  that  the  sign  of  the  slope  is  changed  on 
passing  to  the  M,D-plane.  All  parts  of  Fig.  7 therefore  yield  a curve  of  the  general  shape 
of  that  in  Fig.  6 again,  the  extremum  E corresponding  to  the  maximum.  As  there  the  portion 


beyond  the  crest  presumably  corresponds  to  unstable  deflagration  and  therefore  is  not  part  of  the 

response.  When  there  is  no  extremum  in  the  strip  (6),  the  end  point  lies  beyond  the 

crest  and  there  is  no  acceptable  response  at  all.  We  conclude  that  no  dramatic  change  in 
response  from  that  for  adiabatic  conditions  takes  place,  in  contrast  to  the  pyrolyzing  solid. 

Finally,  we  come  to  the  effect  of  distributed  heat  exchange  on  the  response.  The  analysis 
in  Sec.  6 makes  no  use  of  the  pyrolysis  law,  which  only  serves  to  locate  on  the  response 

curve.  It  is  therefore  equally  valid  for  the  vaporizing  liquid,  the  Claus ius-Clapeyron  law 
serving  only  the  same  purpose.  In  fact  complete  analogy  obtains  when  there  is  no  radiative 

exchange  at  the  surface:  Tg  becomes  a function  of  D alone  and  the  response  may  be  written 

_2 

(22)  M = M„(D)  e_$M  12  where  M.  = T29_1exp(-6/2T  ) 

(J  U 00  « 

since  T is  a constant.  (Otherwise  T and  T are  functions  of  both  M and  D.) 

oo  S 00 

Arguments  similar  to  those  in  Sec.  6 determine  the  shape  of  the  response  as  D -*•  D , 
showing  that  it  depends  on  whether  T^D^)  is  greater  or  less  than  T^.  When  it  is  less 
the  curve  eventually  goes  above  that  in  Fig.  6 and  asymptotes  D = D . The  sudden  complete 
evaporation  is  therefore  replaced  by  an  increasing  rapid  one  as  the  flame  recedes  to  infinity. 
That  was  the  behavior  reported  by  Kapila  and  Ludford  (1977)  who,  by  taking  T T^  ensured 
T (D  ) < T.  because  of  the  endothermic  vaporization.  However,  for  T (D  ) > T.  , the  curve 
continues  to  dip  down  until  the  D-axis  is  approached  when  it  bends  back  to  the  right  (because 
of  the  M in  the  exponential)  to  provide  a C,  exhibiting  true  extinction.  It  does  not 
take  much  imagination  to  foresee  that  an  oval  will  be  formed,  the  part  beyond  the  crest 
closing  with  the  lower  part  of  the  C.  Such  ovals  are  reported  by  Spalding  (i960)  and  are 
sketched  by  Williams  (1965);  they  were  supposedly  valid  for  heat  loss  by  surface  radiation 
alone,  whereas  we  see  that  they  can  only  occur  when  there  is  distributed  heat  loss  (at 
least  when  9 ->•  «>). 

The  picture  near  is  unaffected  by  radiative  exchange  at  the  surface  because, 

locally,  T in  the  formula  (22)  can  be  replaced  by  its  value  at  P . Hence  the  behavior  for 
$ small  is  the  same,  at  least  in  the  neighborhood  of  effective  extinction. 

-22- 
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